Ultra-relativistic nuclear collisions: event shape engineering 
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The evolution of the system created in a high energy nuclear collision is very sensitive to the 
fluctuations in the initial geometry of the system. In this letter we show how one can utilize these 
large fluctuations to select events corresponding to a specific initial shape. Such an "event shape 
engineering" opens many new possibilities in quantitative test of the theory of high energy nuclear 
collisions and understanding the properties of high density hot QCD matter. 
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Many features of multiparticle production in ultra- 
relativistic nuclear collisions reflect the initial collision 
geometry of the system. As the initial conditions affect 
to a different degree all the particles, it leads to truly 
multiparticle effects often referred to as anisotropic col- 
lective flow. Studying anisotropic flow in nuclear colli- 
sions provides unique and invaluable information about 
the evolution of the system created in a collision, prop- 
erties of high density hot QCD matter, and the physics 
of multiparticle production in general [TJ [2]. Recently, 
a significant progress has been reached in understanding 
the role of the fluctuations in the initial density distri- 
bution [3]-[7]. In particular it was realized that such fluc- 
tuations lead to odd harmonic anisotropic flow, which 
enable new insights into dynamics of the system evolu- 
tion. The experimental measurements [51 [S] confirm the 
existence of collective flow up to at least sixth harmonic, 
thus validating the picture. 

At present, the effect of the initial geometry on final 
state observables can be studied only by varying the col- 
lision centrality, or colliding nuclei of different size and 
shape. It has been always tempting to study anisotropic 
flow at maximum particle density, but this was possi- 
ble only in very central collisions where the anisotropies 
are small. Collisions of very non-spherical nuclei, such 
as uranium, should be able to provide events with large 
initial anisotropy and high particle density (in the so- 
called body-body collisions), however the analysis might 
be very complicated due to large variety of possible over- 
lap geometries that have to be experimentally disentan- 
gled. In this paper we discuss how one can select events 
corresponding to different initial system shapes utilizing 
strong fluctuations in the initial geometry even at fixed 
impact parameter, e.g. Au+Au central collisions but of 
large initial anisotropy, and in this way to study the sys- 
tem evolution under conditions not possible before. 

The study of particle production in the events corre- 
sponding to a specific geometry opens a number of very 
attractive possibilities. One of those, mentioned above, is 
the study of the system evolution in a high density regime 
(central collisions) and concurrently strongly anisotropic 
initial conditions. This would add new constraints to 



questions such as the approaching of the system evolution 
toward the so-called "hydrodynamic limit" and the devel- 
opment of the anisotropic flow velocities fields. Analysis 
of transverse momentum spectra in such events can shed 
light on correlation between radial and anisotropic flow. 
Another example would be understanding the "away- 
side" double bump structure in two-particle azimuthal 
correlations [10l Qj]. Several years ago, this attracted 
a considerable attention as a possible indication of the 
Mach cone due to propagation of a very energetic par- 
ton through the dense medium. More recently it was 
found that this structure is likely due to triangular (third 
harmonic) flow. Additional proof for this interpretation 
might come from studying such correlations in events 
with very small triangularity. Several other examples, 
including azimuthally sensitive femtoscopy and an esti- 
mate of the background effects in chiral magnetic effect 
studies will be discussed later in the paper. 

There might be different approaches to perform such 
an event shape engineering (ESE). The one adopted in 
this paper is an an extension of the technique proposed 
in [T2J that is based on the event selection according to 
the magnitude of the so-called reduced flow vector q n 
(the subscript n is the harmonic number, for the exact 
definitions see below) . We always perform ESE using two 
subevents. We use here a common terminology in flow 
analyses, where a subevent refers to a distinct subset of 
all measured particles selected either at random or in a 
given momentum region. One of the subevents is used for 
the event selection (we will always call it subevent "a" 
below) whereas the physical analysis is performed on the 
second subevent (subevent "b"). Using two subevents 
helps to avoid nonphysical biases due to nonflow effects 
as discussed below. We use Monte-Carlo Glauber model 
to illustrate how the event selection based on flow vectors 
works and outline the general scheme for the correspond- 
ing experimental analysis. 

To quantify the anisotropic flow we use a standard 
Fourier decomposition of the azimuthal particle distri- 
bution with respect to the n-th harmonic symmetry 
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planes [HDI]: 
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where v n is the n-th harmonic flow coefficient and ty n 
is the n-th harmonic symmetry plane determined by the 
initial geometry of the system (as given by the partic- 
ipant nucleon distribution, see below). The event-by- 
event fluctuations in anisotropic flow are believed to fol- 
low the fluctuations in the corresponding eccentricities 
of the initial density distribution. Following [7], for the 
latter wc use the definition 

£ n ,x = (r n cos(n0)) , e n>v = (r™ sin(n^)) (2) 

£n,p — \J &n,x ^n,y> tan(?l^l' n ) = £n,y / &n,x ; (3) 

where s n , P is the so-called participant eccentricity. The 
average can be taken with energy or entropy density as 
a weight. In our Monte-Carlo model we weight with the 
number of participating (undergoing inelastic collision) 
nucleons. For the nucleon distribution in the nuclei we 
use Woods-Saxon density distribution with the standard 
parameters (for the exact values see [H]); the inelastic 
nucleon-nucleon cross section is taken to be 64 mb. We 
assume that the flow values are proportional to the cor- 
responding eccentricities with the ratio fixed to approxi- 
mately reproduce measured v n values [9 . As it is shown 
in |15) . in this case the distribution in v n is very well 
described by the so-called Bessel-Gaussian (BG) distri- 
bution BG(w; «o, (Jvx), where 
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which is a radial projection of 2-dimensional Gaussian 
distribution with the width a in each dimension and 
shifted off the origin by distance x$. 

The flow vectors are calculated in two subevents [TJ Q3] 
with multiplicities in each subevent approximately cor- 
responding to An = 0.8 in Pb+Pb collisions at LHC 
energies [16] (approximately 1200 charged particles per 
subevent for 0-5% centrality). The multiplicities are gen- 
erated with a negative binomial distribution based on 
number of participants and number of the binary colli- 
sion as in [13] . For each 5% width centrality bin discussed 
below, we analyze about 1.2 M simulated events. 

The flow vectors are defined as 
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Qn,x = y^cos(n<ftj); Q, hy = y^sin(n^); 



q n = Qn/vM, 

«£ = l + (M-l)<coB[n(&-^)]> 
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where M is the particle multiplicity and <pi are the par- 
ticle azimuthal angles of particles in a given subevent. 
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FIG. 1. (color online) Mean elliptic and triangular flow values 
in a-subevent as function of the corresponding q n magnitude 
in 6-subevent. 



Eq. [7] presents the relation of the length of q n vector to 
the average correlation between all pairs of particles in a 
given event. The event-by-event distribution in the mag- 
nitude of flow vectors q„ has been proposed 13 and often 
used to measure the average flow [1] [17] . The distribu- 
tion in q n is determined by the v n distribution convoluted 
with statistical fluctuations due to finite multiplicity. For 
relatively high multiplicities (M > 300) it is very well 
described by BG distribution BG(g; qo, a qx ) with param- 
eters related to those of v„ distribution: 
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= VMv , a% = - [1 + (M - l)(2c4 + 5)] , (8) 



where M is the multiplicity used to build the flow vector, 
and a nonflow parameter 8 accounts for possible correla- 
tions not related to the initial geometry of the system. 
(For a more detailed discussion of the functional form of 
q n distributions see [18].) Thus, the fit to ^-distribution 
provides information about underlying flow fluctuations, 
if the nonflow contribution can be neglected or estimated 
from other measurements. 

Zero nonflow. We start the discussion of the ESE with 
the simplest case when all the correlations in the sys- 
tem are determined only by anisotropic flow. Figure [l] 
shows the average values of calculated via 2-particle 
correlation method in one of the subevent ( "b" ) as func- 
tion of the flow vector magnitude in the second subevent 
( "a" ) . We remind the reader, that in this simulations the 
two subevents are statistically independent and are cor- 
related only via common participant plane and flow val- 
ues. There are no nonflow correlations included at this 
stage. In this case the results for b {2} coincide with 
"true" values of (vfy (not shown), though have slightly 
larger statistical errors due to finite multiplicity of the 
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subevent. 

The results in Fig. [T] demonstrate, that depending on 
<?2, a one can select events with average flow values varying 
more than a factor of two. How well one can "resolve" the 
flow fluctuations depends on the number of particles used 
to calculate the flow vector as well as, though weakly, on 
flow magnitude itself. We find that for centrality 20-25%, 
the width of the vi distribution for a fixed (72, a value is 
about factor of 1.5 smaller than that for unbiased event 
sample (changing from 0.031 to 0.022); it decreases for 
about 20% if one double the size of the subevent (double 
the multiplicity) used for 92 determination. 

Let us demostrate now how in practice one can obtain 
an information about the v n distributions, correspond- 
ing to different cuts on the q na values, from the fits to 
the g^b-distributions. Figure [2] shows distribution in q n ^ 
(subevent-b) for three different cuts on q n ,<n separately 
for the second and third harmonic flow. All q n t, distribu- 
tions in Fig.[2]are fit to the BG functional form to extract 
the corresponding mean flow values and the correspond- 
ing width (see, e.g. [1]). It is remarkable that the fits 
are very good not only for the unbiased g-distributions 
but also to the ones corresponding to the low flow and 
high flow "engineered events" (corresponding to the 5% 
lowest and 5% highest q n , a events). Using the extracted 
parameters we plot the corresponding v n distributions in 
Fig. [3] (shown by dashed lines) and compare to the actual 
("true") v n distributions, which is known in this Monte- 
Carlo simulation (shown as a histogram). One finds an 
excellent agreement between the two indicating that the 
v n distributions in the "shape engineered" events are very 
close to the BG form. 

Nonflow effects. The ESE approach described above is 
based on using two subevents. In this case possible non- 
flow effects can be separated in two major categories 
(a) when nonflow effects are present within each of the 
subevents, but there is no nonflow correlations between 
subevents "a" and "b", and (b) when nonflow correla- 
tions are present between, as well as within, subevents. 
As we show below one should try to minimize the nonflow 
correlations between the two subevents which are used for 
ESE selection and physical analisys, respectively. A prac- 
tical solution to that might be to use subevents which are 
separated by a significant (pseudo)rapidity gap. 

The case (a) does present a certain challenges to the 
analysis, but no more than the one in the conventional 
flow analysis. Once the event selection is done with q n ^ a 
cuts, the flow in the selected events can be estimated 
using particles in subevent "b" with standard methods 
including many-particle cumulant analysis. The case (b) 
is significantly more complicated. Below we only discuss 
possible biases, without trying to resolve the problem. 

We simulate nonflow effect by assuming that half of all 
particles in the entire event are produced in pairs with 
both particle in a pair emitted with the same azimuthal 
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FIG. 2. (color online) 52,6 and q^^ distributions in the event 
samples selected by different cuts on the corresponding q n , a - 
vector magnitude indicated in the plot. The lines show the 
BG fit to the distribution. 



angle. Each particle is assigned randomly to one of the 
two subevents. In this case the nonflow parameter 6 = 
1/(2M), where M is the (full) event multiplicity, which 
roughly corresponds to the nonflow estimates in real LHC 
events for particles at midrapidity \r]\ < 0.8. 

Figure [4] presents the results for flow calculation in 
subevent "a" using 2- and 4-particle cumulant methods 
as function of q2,b- The expectations based on simulated 
flow are also shown. One observes a significant bias due 
to nonflow, leading to overestimate the flow values in high 
flow selected events and underestimate in the low flow se- 
lected events. This trend is due to positive character of 
the nonflow correlations. The corresponding bias in cor- 
responding v distributions is shown in Fig. [5] Note that 
even though the bias for mean values of flow is somewhat 
modest, at large values of v n the actual distribution could 
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FIG. 3. (color online) Actual (true) V2 and V3 distributions 
in the event samples selected by different cuts on the corre- 
sponding q n -vector magnitude indicated in the plot compared 
to that extracted from the BG fits to q n ,b distribution shown 
in Fig. [2] (dashed lines). Note that the lines are not the fit to 
the histograms! 

differ by order of magnitude from the one deduced from 
(/-distribution fits. 

Below we discuss very briefly several analyses, which 
can profit from the event shape engineering. 

The chiral magnetic effect proposed in |19H21) is a charge 
separation along the magnetic field. A correlator sensi- 
tive to the CME was proposed in Ref. [22] : 

(cos(0 Q + <f> p - 2f RP )) , (9) 

where subscripts a, (3 denotes the particle type. The 
STAR [23J [24], as well as the ALICE [25] collaboration 
measurements of this correlator are consistent with the 
expectation for the CME and can be considered as evi- 
dence of the local strong parity violation. The ambiguity 
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FIG. 4. (color online) Elliptic flow measured with 2- (red 
points) and 4-particle (blue) cumulant method in subevent 
"a" as a function of the corresponding q2,b magnitude. Solid 
symbols correspond to centrality 20-25%, and open symbols 
to 0-5% centrality. The true (simulated) values are shown 
by green markers, as expected for 2-particle cumulant results 
and by magenta for 4-particle cumulant results. 



in the interpretation of experimental results comes from 
a possible background of (the reaction plane dependent) 
correlations not related to CME. Note that a key ingre- 
dient to CME is the strong magnetic field, while all the 
background effects originate in the elliptic flow [22 . This 
can be used for a possible experimental resolution of the 
question. One possibility is to study the effect in central 
collisions of non-spherical uranium nuclei [12] , where the 
relative contributions of the background (proportional to 
the elliptic flow) and the CME (proportional to the mag- 
netic field), should be very different in the tip-tip and 
body-body type collisions. The second possibility would 
be to exploit the large flow fluctuations in heavy-ion colli- 
sions as discussed in [T2J[5S] and the ESE would be a tech- 
nique to perform such an analysis. (Note also that the 
magnetic field depends very weakly on the initial shape 
geometry fluctuations [35].) Yet another test, proposed 
in [27], is based on the idea that the CME, the charge sep- 
aration along the magnetic field, should be zero if mea- 
sured with respect to the 4-th harmonic event planes, 
while the background effects due to flow should be still 
present, albeit smaller in magnitude (~ V4). An exam- 
ple of such a correlator, would be (cos(2</> Q + 2<frp — 4^4), 
where ^4 is the fourth harmonic event plane. The value 
of the background due to flow could be estimated by 
rescaling the correlator Eq. [9j Such measurements will 
require good statistics, and strong fourth harmonic flow. 
Again, the ESE can be very helpful to vary the effects 
related to flow. 



5 




0.05 0.1 0.15 0.2 0.25 0.3 

v 3 

FIG. 5. (color online) The same as in Fig [3] but for the case 
of nonflow described in text. 

Measuring the shape and freeze-out velocity profile with 
azimuthally sensitive femtoscopy. Different shapes in the 
initial geometry of the collision, to a different degree will 
be preserved in the system freeze-out shapes. It was 
shown in [25] that those shapes can be addressed exper- 
imentally with azimuthally sensitive femtoscopic analy- 
sis (HUSO], which has a goal to obtain the geometry of the 
source relative to different harmonic symmetry planes. 
Such an analysis would definitely profit from event with 
extreme values of anisotropy provided by the ESE, as the 
variation of femtoscopic parameters with azimuth would 
be better pronounced. General details of femtoscopic 
analyses and discussion of the experimental results can 
be found in a review [3Tj . 

Summary. Event shape engineering, providing possibility 
to study events corresponding to nuclear collisions with 
different initial geometry configuration, promises wide 
use in further studies of the properties of the strongly 



interacting matter. 
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